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Studies of Fermi data indicate an excess of GeV gamma rays around the Galactic center (GG), 
possibly due to dark matter. We show that young gamma-ray pulsars can yield a similar signal. First, 
a high concentration of GG supernovae naturally leads to a population of kicked pulsars symmetric 
about the GG. Second, while very-young pulsars with soft spectra reside near the Galactic plane, 
pulsars with spectra that have hardened with age accumulate at larger angles. This combination, 
including unresolved foreground pulsars, traces the morphology and spectrum of the Excess. 


Introduction. — The Fermi Large Area Telescope [T] 
has transformed GeV gamma-ray astrophysics. Pulsar 
physics in particular has experienced an enormous ad¬ 
vance [2]. The impact is much greater than even the 
impressive increase in the number of gamma-ray pulsars, 
from 7 pie-Fermi to >170 now HE], alone suggests. GeV 
detections of old, recycled millisecond pulsars (MSPs) 
went from zero to ^ 60 in under six years. The >100 
detections of young pulsars (t < Myr) imply ubiquitous 
GeV emission EllS], with the radio-quiet population now 
at ~ 40 from a lone prototypical member, Geminga mM- 

Fermi could also fulfill the long-standing hope of de¬ 
tecting gamma rays from the annihilation or decay of 
dark matter im, which comprises « 84% of all mat¬ 
ter [HUH]. A number of groups [HHZil have used Fermi 
data to search for a signal originating from around the 
Galactic center (GC), where this flux should be largest. 

The exciting recent development from these studies is 
the measurement of an extended excess of GeV gamma 
rays above model predictions peaking from 1-3 GeV. The 
Excess spectrum is reasonably fit by the annihilation of a 
30—60 GeV dark matter particle [e.g.. [HHH] . though the 
angular intensity implies a density profile steeper than 
often extrapolated from simulations. 

The great challenge is to determine whether this sig¬ 
nal indeed arises from dark matter or is due to an 
unaccounted-for source of gamma rays. One suggestion 
is a contribution from MSPs [HI SSHIHI : which have a 
similar spectral shape as the Excess in gamma rays and 
are expected to exist there at some level [49] . It has been 
argued that MSPs cannot match the Excess [snun] and 
uncertainties in MSP formation make definite predictions 
difficult. 

We show that young pulsars from the GC itself produce 
a diffuse GeV flux that has been widely underestimated. 
Such pulsars arise from core-collapse supernova explo¬ 
sions of short-lived massive stars. Their birthplaces thus 
tend to trace star formation. The GC is the most con¬ 
centrated star forming region in the Milky Way (MW); 
the inner ^ 200 pc central molecular zone (CMZ) ac¬ 


counts for 5 — 10% of the current Galactic star formation 
rate [e.g., IKHUHH] . Further, the CMZ contains an esti¬ 
mated ~ 13% of all MW Wolf rayet stars [57U6T] : these 
evolved, 25 Mq stars are near explosion and may re¬ 
main from an even more intense recent period of star 
formation [e.g., [52] ■ This implies a substantial, ongoing 
production of pulsars, with birth kicks leading to a con¬ 
tinuous, symmetrical distribution centered on the GC. 

Fig.g displays a second piece to this puzzle. Here 
we compare the gamma-ray spectrum measured from 
Geminga, the best characterized of radio-quiet pulsars, 
to that expected from 35 GeV dark matter annihilating 
to a 6-quark pair [63] as proposed for the Excess. The 
striking similarity to the Excess is a consequence of the 
hardening of gamma-ray pulsar spectra with decreasing 
spin-down power (E) as seen in Fermi data [2| El] . 



FIG. 1: A comparison of phase-averaged Fermi gamma-ray 
data of Geminga |52| to our spectral shape from pulsars within 
5° of the GC and > 2° from the disk (see text), which resemble 
that from 35 GeV dark matter particles annihilating to bb as 
proposed to explain the GC Excess [25]. 
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Very young, high-i? pulsars with soft spectra then re¬ 
side within the Galactic plane, where they are masked 
or subtracted as point sources and do not contribute to 
the observed Excess. On the other hand, pulsars old 
enough to reach large angles have hard spectra (Geminga 
has a characteristic age of ~3xl0®yr). Fermi has also 
shown that the gamma-ray luminosity, oc de¬ 

clines much more slowly than spin-down power, while 
gamma-ray efficiency, L-^/E, rises at low E [3]. Since an 
additional 20% of Galactic star formation occurs within 
a projected radius of 15° of the GG [5S], there should be 
unresolved high-latitude disk pulsars that also contribute 
to the Excess. Indeed, ~ 15% of all radio pulsars are seen 
in this region [66]. 

Our conservative estimates suggest the population of 
unresolved young pulsars may constitute a significant 
fraction, if not all, of the excess GeV gamma rays. 
Thanks to intensive Fermi and radio studies, we are able 
to forward model the emission of young pulsars. This is in 
contrast to modeling MSPs, which have complicated for¬ 
mation channels, or ancient cosmic-ray bursts [671168j or 
dark matter which have a large range of free parameters. 
We do not include any MSP or other non-young-pulsar 
contribution in this work. We discuss a few of the many 
implications here and in a companion paper |69j . 

The GC pulsar factory. — We use a suite of Monte 
Garlo simulations to forward model the present day dis¬ 
tribution of ordinary young pulsars throughout the GG 
and Galactic disk. Our methods largely follow Refs. HQ] 
and [5], who successfully reproduced the observed distri¬ 
bution of nearby radio and gamma-ray pulsars, respec¬ 
tively. We extend both by including the GMZ contribu¬ 
tion (avoiding here the peculiar central parsec [7T]~1. 

To generate our simulated gamma-ray pulsar popula¬ 
tion, we select the pulsars’ ages randomly assuming a 
constant birth rate of 2.1 x 10“^ yr“^ in the disk HQ] [72], 
formed on circular orbits within the MW spiral arms m- 

The GG rate depends on the number of massive stars 
formed over the previous ~ 30 Myr, due to their finite 
ages m, while the present star formation rate may be 
near a minimum if episodic on ~ 10 Myr periods |62j . 
Over this time their positions will be symmetrized by 
cluster dispersion and differential rotation. For simplic¬ 
ity, we assume a constant birth rate of 1.5 x 10“^yr“^ 
in the GMZ, « 7% of the Galactic rate. These pul¬ 
sars are formed from a spherical isothermal distribution 
(n oc r~^), between 20 and 200 pc of Sgr A* - the region 
of most GG star formation. 

We kick each pulsar by selecting a random kick veloc¬ 
ity from a three-dimensional normal distribution with a 
mean velocity of « 408kms~^ [74]. Each pulsar orbit is 
integrated over its lifetime using GALPY [106] with the 
MWPotential2014 model for the Galactic potential m- 

The pulsar gamma-ray luminosity depends on its total 
spin-down power, which in our models is determined by 
the surface magnetic field, initial spin period, and age. 



b (degrees) 


FIG. 2: Fluxes of simulated pulsars within Galactic longitude 
|1| < 20° and latitude |6| < 20° from one realization of our 
fiducial model, separated into those from the GC (black di¬ 
amonds) and disk (small red dots). These are compared to 
the population of Fermi 3FGL point sources (blue circles), 
gamma-ray discovered pulsars (red stars), and radio discov¬ 
ered pulsars (green triangles) [3l|4l[76]. The dashed line shows 
the 50% completeness limit for 2FGL [see our Eq. ]2f;|Z7]. 

We assign each pulsar a constant surface B from a log¬ 
normal distribution with mean (logj^Q B) = 12.65 and 
<^\ogB = 0.45. The variance in the magnetic fields is 
roughly between those used by Refs. [7D] & [5]. The 
pulsars are born with an initial spin period, Pq following 
[^. Therefore the spin-down luminosity of the pulsar 
E = 47r^JPP“^, where 7 = 1 x lO^^gcm^ is the pulsar 
moment of inertia, P = 87r^7?®i?^/(3/c^P), R = 10® cm 
is the radius of the pulsar, c is the speed of light, and 
P = (p2 + 167r2p®p2^/(3/c®))i/2. 

Gamma rays in excess.— Fermi observations of 
young pulsars show that their sky-averaged gamma-ray 
luminosity follows a trend that agrees with the theoreti¬ 
cally expected L.y oc [e.g., [21 Oj although with no¬ 

table scatter. We normalize for each pulsar as 

= (is4i7;) w 

The heuristic value of C = 1.0 mm underestimates the 
luminosity of Geminga by a factor of 5, and the mean and 
median gamma-ray luminosity of ordinary pulsars found 
in [3] by a factor of « 3. For a given pulsar, models pre¬ 
dict different beaming corrections to the observed flux 
dependent upon the angle to the observer and magnetic 
field alignment m- Using the beaming-corrected rela¬ 
tions of [64] modeled from many individual light curves, 
we find that the mean luminosity is C = 2.2, after re¬ 
moving the top and bottom 5% of the sample. Here, we 
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FIG. 3: The spatial distribution of 2 GeV gamma rays from 
unresolved young pulsars in one realization. Note that we do 
not include detector noise or other Galactic emission here. 


are mostly concerned with a population-averaged lumi¬ 
nosity density and take a fiducial value of C = 1.3, which 
produces one-third of the unassociated plane sources in 
the Fermi 3FGL catalog m- For simplicity, we assume 
that gamma-ray emission turns oS at E ^ erg s“^, 

higher than the outer gap death line of [sniisi]. 

In Fig. we show the gamma-ray fluxes of a single 
realization of the GC and Galactic disk pulsar popula¬ 
tions as a function of Galactic latitude, b. Also plotted 
are the fluxes of 3FGL point sources m and gamma-ray 
pulsars [311 Eg within the same region. In our model, 
no pulsars formed at the GC are bright enough to be de¬ 
tected as a distinct point source by Fermi. Indeed, only 
pulsars from the Galactic disk cross the 50% point source 
completeness limit we later use for subtraction. 

This same region contains seven gamma-ray-discovered 
pulsars and our simulations typically contain a compara¬ 
ble number. We also find that the results of our sim¬ 
ulations are compatible with the entire Galactic pulsar 
population. Like previous analyses Oil], we produce too 
few of the most luminous pulsars. Increasing C helps to 
alleviate some of this tension (see [69] for a discussion 
of this and pulsar beaming). Fig. [2 also shows that, on 
average, the highest GC pulsars are the closest to the 
GC, a natural consequence of spin-down. In our model, 
pulsar gamma-ray luminosity decreases with time as . 

To generate gamma-ray maps of the GC, we convolve 
each pulsar with the Fermi point spread function (PSF) 
[see Eq. (5) of Ref. |H3 twice |107j . We use luminosity- 
weighted spectra of the pulsars to estimate their flux den¬ 
sity at 2 GeV, as detailed in [69]. In this analysis we re¬ 
move all pulsars with ages t < 10"^ yr, as the gamma-ray 
luminosity of the youngest pulsars is poorly understood. 
To generate the diffuse map, we then subtract sources 


with fluxes greater than the 50% completeness limit of 
Fermi derived from the brightness distribution of 2FGL 
point sources between |5| < 20° and |Z| < 20° [77] 

Lp, = 1 X ergs-i cm'^, (2) 

where b is the Galactic latitude in degrees, as shown in 
Fig. We use 2FGL only to be consistent with p51[?7] . 

In Fig. [^ we show the surface brightness map of a sin¬ 
gle realization of the unresolved gamma-ray pulsars to¬ 
wards the GC. The flux is centrally concentrated around 
the CMZ, but outside of the central few degrees, pulsars 
in the Galactic disk become the dominant component, 
with a population of point sources extending well above 
the Galactic plane. The outermost regions are the most 
impacted by Poisson fluctuations in the underlying pulsar 
population, and show the greatest amount of asphericity. 
Since most of the flux in the inner 5° comes from the 
CMZ pulsars, flux-weighted measures of the asymmetry 
of the emission will not be heavily impacted by the disk. 

In Fig. [^ we compare the surface brightness profile of 
young pulsars to the gamma-ray excess at 2 GeV reported 
by Daylan et al. [25] and Galore et al. [23 [28]. The solid 
red line shows the mean surface brightness of the pulsars 
using concentric rings with 1° width that exclude the 
Galactic plane (at |6| < 2°), motivated by comparison 
with [23, while the green dashed line includes only the 
high-latitude pulsars (|5| > |Z| and |6| > 2°). 

We see that the total flux of the young pulsars can 
plausibly explain much of the GeV Excess in the region 
we explore in this work, with a slope well within the liter¬ 
ature range. Since pulsar behavior near the gamma-ray 
death line remains unclear, as is the underlying mag¬ 
netic field distribution of the pulsars, we have chosen to 
plot along with the total flux implied by Eq. 0 down 
to if > erg s“^ (solid red line) the range of un¬ 

certainty of the underlying model (gray band). We have 
done this by varying the main parameters within the em¬ 
pirically observed limits. We have varied the cutoff be¬ 
tween 1034 > E > 1033 erg s-^, 0.55 > uiogs > 0.30, 
and 2.2 > (7 > 1.0. Currently, the lowest E for a young 
gamma-ray pulsar is ~ 3x 1033 erg s”^ [g. Unfortunately, 
low-if radio-quiet pulsars often lack distances (thus lumi¬ 
nosities). Indeed, tests coupling both local observations 
and the large-GC angle pulsar population may improve 
our understanding of these faintest pulsars [3E61IH3. 

The inferred spectrum of the Excess mostly comes from 
the ~ 5° region around the GC [23 [2H]- In our model, the 
luminosity from this region is dominated by Geminga-like 
pulsars (which has E « 1034-3 erg s“4). Fig.j^shows the 
flux-weighted spectrum of all pulsars that are beyond the 
plane (|6| > 2°) and within 5° of the GC, which is com¬ 
parable to the Excess spectral shape near the peak. Note 
that ^ 10 GeV pulsar emission may be underestimated 
by using an exponential cut-off [e.g., |3 US IM]- 
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FIG. 4: Intensity of 2 GeV backgrounds as a function of angular separation from the Galactic Genter. The solid red line shows 
the average pulsar contribution from our fiducial model, obtained from angular rings (excluding the |&| < 2° disk), with the 
gray shaded region showing the systematic uncertainty in the gamma-ray pulsar population. The green dashed line is only from 
pulsars at high Galactic latitude (|fe| > |1| and |6| > 2°). Also shown are intensities of the gamma-ray Excess, from Ref. |25| 
which used Gaussian angular rings (masking |&| < 2°) out to « 10° (black points and error bars; the lower point at 5° is their 
“Full Sky”, similar to their “Galactic Genter” values) and Ref. |28| (shaded blue bands) which assume a fixed profile shape. 


Discussion. — We have shown that a population of 
young ordinary pulsars originating in the central molec¬ 
ular zone of the Galactic center and the Galactic disk 
naturally produces a gamma-ray spectrum, amplitude, 
and spatial distribution resembling the observed GeV ex¬ 
cess in the Galactic center. Our model uses the observed 
properties of radio [701 m] and gamma-ray pulsars |5] 
to estimate their gamma-ray flux. The young pulsar 
luminosity profile is then determined by kicked pulsars 
evolving with age, while we include no contribution from 
millisecond pulsars. 

Galore et al. [75] found that the Excess brightness var¬ 
ied between different analyses at the 3—tr level, suggest¬ 
ing at least a factor of two uncertainty in its absolute nor¬ 
malization. In this work, we used a common, simplified 
model for pulsar evolution. This model also has a factor 
of two uncertainty due to discrepancies in the preferred 
surface magnetic fields for radio pulsars and gamma-ray 
pulsars, pulsar behavior near the gamma-ray death line, 
and the GC SN rate (though our fiducial value is near 
that preferred in Fermi Bubble models [511[5S]). We dis¬ 
cuss these in [55] , although we note here that the param¬ 
eters used by our models were specifically optimized to 
reproduce the currently observed distribution of pulsars. 

Nevertheless, there remains much room for improve¬ 
ment. While current “gap” models of pulsar emission 
[e.g.. ISBHOO] work reasonably well [64jj ab initio models 
under development promise better physical understand¬ 
ing [5TH51] . Ideally, we would also have a similar under¬ 
standing of neutron star birth kicks, spins, masses, and 
magnetic fields [55H55] . 


We discuss further observations to help determine the 
young pulsar fraction in |69j . These include the pul¬ 
sar contribution to a similar GeV excess seen along the 
Galactic disk [55], which may also be evident as residual 
emission in the inner Galaxy |29] . We have not here in¬ 
cluded electrons and positrons from pulsar winds (e.g., 
Geminga displays diffuse TeV emission [lOOj due to 
[lOlj ). If proportional to E, most of the output is con¬ 
centrated near the plane to yield an inverse Compton flux 
(as may be needed in the inner galaxy mm)- 

Despite present uncertainties, young pulsars must be 
accounted for in future models of the Excess, with our 
conservative estimates suggesting that they constitute a 
substantial fraction. If the Excess extends to higher en¬ 
ergies, as in [75], and if pulsars only contribute a fraction 
at GeV energies, a higher-mass dark matter candidate 
may yet be present (see [55]). Unlike the GC, old dwarf 
galaxies should not contain young pulsars, so improved 
gamma-ray studies [14] I102j will help our understanding 
of both the Excess and pulsar physics. Clearly, additional 
Fermi data will be vital to resolve more sources, charac¬ 
terize pulsars, and accumulate statistics at the highest 
energies, where upcoming lACTs will also be relevant 
[103H105] . to dissect the Excess. 
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